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Abstract 



We provide a new proof for the global well-posedness of systems coupling 
fluids and polymers in two space dimensions. Compared to the well-known ex- 
isting method based on the losing a priori estimates, our method is more direct 
and much simpler. The co-rotational FENE dumbbell model and the coupling 
r^ , Smoluchowski and Navier-Stokes equations are studied as examples to illustrate 

Cd ' our main ideas. 
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c^j . 1 Introduction 

en 

^ I The dynamics of many polymeric fluids are described by two-scale micro-macro mod- 

els. The systems usually consist of a macroscopic momentum equation and a micro- 
scopic Fokker-Planck type equation. The fluid is described by the macroscopic equation 
(sometimes the incompressible Navier-Stokes equations), with an induced elastic stress. 
The stress is the micro-macro interaction. The particles in the system are represented 
by a probability distribution ip(t,x,R) or ip(t,x,m) that depends on time t, macro- 
scopic physical location x and particle configuration R or m. The Lagrangian transport 
of the particles is modeled using a Taylor expansion of the velocity field, which accounts 
for a drift term that depends on the spatial gradient of velocity. The system attempts 
to describe the behavior of this complex mixture of polymers and fluids. For more 
physical and mechanical backgrounds, see [H ITT] . 
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The FENE (Finite Extensible Nonlinear Elastic) dumbbell model is one of the typ- 
ical and extensively studied micro-macro models. In this model a polymer is idealized 
as an "elastic dumbbell" consisting of two "beads" joined by a spring which can be 
modeled by a vector R. Mathematically, this system reads 

dtu + {u ■ V)u + Vp = uAu + div r, x e M^, 

div M = 0, X e R^, 

dt^ + {u ■ V)^ = divR [ - W{u) ■ R^ + /3Vr^ 
+VRUij] , (x, R)eR'^ X 5(0, Ro), 



:i.ii 



{VrU^ + /3Vr^) ■ n 
t = : u(t,x) = Uo{x) 



0, ondB{0,Ro), 

ipit, X, R) = ipo{x, R)- 



In the above system, u = u(t,x) denotes the velocity field of the fluid, p = p(t,x) de- 
notes the pressure, iplt, x, R) is the distribution function for the internal configuration, 
i^ > is the viscosity of the fiuid and /3 is related to the temperature of the system. 
Moreover, the spring potential U and the induced elastic stress r is given by 



U{R) = -A;ln(l-|i?| 



|i?nP) 



Ti, 



{R,®VRU)ij{t,x,R)dR. (1.2) 



B(0,i?o) 



Here A; > is a constant. The boundary condition insures the conservation of the 
polymer density. Assume that W{u) = "~^ "^ , which corresponds to the co- rotational 
case. For simplicity of writing, assume that (3 = 1 and Ro = 1 and denote B{0,Ro) by 
B. In what follows, without special claim, V represents V^^, and div represents div^;. 

The Smoluchowski equation coupled with the incompressible Navier-Stokes equa- 
tions is another extensively studied micro-macro model. The Smoluchowski equation 
describes the temporal evolution of the probability distribution function ip for direc- 
tions of rod-like particles in a suspension. Mathematically, the system reads 

' dtU + {u ■ V)u + Vp = uAu + div r, x e M^, 

div M = 0, X e M^, 

dttp + (u ■ V)^ + diVg{G{u, ^)V') - Ag^ = 0, (x, m) eR'^ x M, 
t = : u(t,x) = Uo{x), ipit^x^m) = tpQ^x^m). 



:i.3) 



Here M is a d-dimensional smooth compact Riemannian manifold without boundary of 
and dm is the Riemannian volume element of M, u = u{t, x) and p = p{t, x) denote the 
velocity field and the pressure of fiuid respectively, ip = tpit^x^m) is the distribution 
function, G{u,iIj) = VgU + W stands for a meanfield potential resulting from the 
excluded volume effects due to steric forces between molecules with W = d^{m)djUi. 



Besides, the added stress tensor r and the potential U are given by 

U{t,x,m)= / K{7n,q)'ijj{t,x,q)dq, 
Jm 

Tij(t,x)= / 'ylj{m)ip{t,x,m)dm (1.4) 



M 

% 
M JM 

Here the kernel if is a smooth symmetric function defined on M x M. 7|- and 7}- 
are smooth, time independent and x independent. 

At present there have been extensive and systematical studies on the existence and 
regularity theories of those 2D micro-macro models of polymeric fiuids [21 El H El [101 
1201 [23]. For example, the first global well-posedness result for FENE (II. ip was derived 
by Lin, Zhang and Zhang [22j for k > 6. Masmoudi [23] extends it to the case of fc > 
by a crucial observation on the linear operator. Very recently the global existence of 
weak solutions to the FENE dumbbell model of polymeric fiows for a very general 
class of potentials was also obtained by Masmoudi [21] . The global well-posedness of 
nonlinear Fokker-Planck system coupled with Navier-Stokes equations (11.31) in 2D has 
been proven by Constantin and Masmoudi in [TO]. When the nonlinear Fokker-Planck 
equation is driven by a time averaged Navier-Stokes system in 2D, global well-posedness 
has been obtained by Constantin- Fefferman-Titi-Zarnescu [9]. 

Most proofs of the above global well-posedness theorems are based on an important 
analytic technique called losing a priori estimate in the spirit of Bahouri-Chemin [2j 
and Chemin-Masmoudi [7|. In [T7j, we studied the blow-up criteria of a macroscopic 
viscoelastic Oldroyd-B system avoiding using the losing a priori estimates. The main 
purpose of this paper is to extend the method in [17] to micro-macro models and 
provide a new proof for global well-posedness of the co-rotational FENE dumbbell 
model (11.11) and the coupling Smoluchowski and incompressible Navier-Stokes equations 
(ll.3p . Compared to the proofs of the theorems in [23] and [TO] , which are based on the 
technique of losing a priori estimates, ours are direct and much simpler. 

For the co- rotational FENE model (II. ID . we have 

Theorem 1.1. Assume that uq G H'{R'^),{s > 2), tpo G H'{R'^;C' f] C^) for some 
r > 2 such that (r — l)k > 1 with J^tpodR = l,a.e. in x. Then there exists a unique 
global solution {u, ^p) of the FENE model ^^D m C([0, 00); iJ^) xC([0, cx)); E'i^^; C?)). 
Moreover, u G ^^,^(0, 00; E'^^) and ^ G Lf^^iQ, 00; H'iR?; C^'')). 

For the definition of £*" and C^'^, please refer to section 2. 

Similarly, for the coupling Smoluchowski and Navier-Stokes system (II. 3p . we have 

Theorem 1.2. Take Uq G W^+''''{R'^) n ^^(M^) and ^0 e W^'''{R^, H-'{M)), for some 
r > 2, e > 0, s > ^ + 1 and ^po > 0. ^j^ilJodm G L\R?) n L'^iR?). Then 1113) has a 
global solution mu e L^,(0, 00; W^'')nLl^{0, cx); W^'"^) andtp G L~ (0, 00; W^'^'iR?; H~'{M))). 
Moreover, for T >To>0, we have u G L~((To,T); W^'"'''). 



We end this introduction by mentioning some other results on micro-macro models. 
Global existence of weak solutions can be found in |U |2T] and local existence of strong 
solutions are studied in [121 [ISl 125]. For macroscopic models, we refer the readers to 
[H [la [Tg US [m E] as references. 

The paper is organized as follows. In section 2, we give some preliminaries. Then 
we give some a priori estimates for the FENE model (11. ip in section 3 and the coupling 
Smoluchowski and Navier-Stokes equations (11. 3p in section 4. The a priori estimates 
obtained in sections 3 and 4 are enough to get the global existence of systems (II. ip 
and (II. 3p and prove Theorem 1.1 and Theorem 1.2 p3| [TO]. 

2 Definitions and Useful Lemmas 

We will use the Littlewood-Paley decomposition in the following sections. Define C to 
be the ring 

C = {ie^':\<\i\< ^}, 
and define V to be the ball 

Let X s-^d if be two smooth nonnegative radial functions supported respectively in V 
and C, such that 

X(0 + Y. ^'(^"'O = 1 for e e M^ and J^ ip{2-'^i) = 1 for ^ G M' \ {0}. 

g>0 



Let us denote by J-" the Fourier transform on M? and denote 

h = T~^(p, h = T^^x- 
The frequency localization operator is defined by 

A,M = T-^ {^{2-'^i)T{u)\ = 2^'? / hi2'^y)uix - y)dy, 

and 

S,u = J'-' [x{2-'0J'iu)] = 2'' [ h{2^y)u{x - y)dy. 

Hence, for s < 2/p or s = 2/p and r = 1, the homogeneous Besov space 5*^ is 
defined as the closure of compactly supported smooth functions under the norm 11 ■ 11 os , 



p,r 



B^r ~ ll*^-^ \\^<l"'\\LP)qeZ\\lr( 



when p = r = oo, s = k + a where k E N and a G (0, 1), then 5^^ turns to be the 
homogeneous Holder space C*^"*"". 



Another kind of space to be used is L^iti, t2] C""), which is the space of distributions 
u such that 



ll^llLP(ti,t2;C"-) 

For the FENE model, let 



S\lY>2'^''\\/\qU\\Lv(ti,t2;L^) < OO- 



-W(fl) 



/^ e-"(«)rf/? 



cd- \R 



2\k 



where the constant c is given such that Jg tp^odR = 1. In fact, {u, ip) = (0, ipoo) defines 
a stationary solution of fll.ll) . 

For r > 1, denote C^ and C^'^ the weighted spaces 



C 



i,: 



£'■ 



^o 



B 



^ 



'l,r 



^ G £^' : 



V^c 



B 






di? < OO 



^c 



di? < OO 



We will need to use the following well-known inequalities. 

Lemma 2.1 (Bernstein Inequalities |5j). For s &M., 1 < p < r < oo and g G Z, one 

has 



l^gWlU' 



<C-2''^p-r^''\\Agu\\LP 



C2'3'\\AgU\\LP < \\\/'AgU\\LP < C2'"\\AgU\\LP, 

\\\V\'SgU\\Lr<C2'''\\u\\Lr, 

ce-^^""\\Agu\\Loo < \\e'^Agu\\Loo < Ce-'^^'^WAguHLo.. 
Here C and c are positive constants independent of s,p and q. 

We will also need the following lemma, whose proof can be found in [T7] . 

Lemma 2.2. Assume that (3 > 0. Then there exists a positive constant C > such 
that 

\\Vg{s,-)\\L^ds<c(l + j \\g{s,-)\\L2ds 

+ sup/ \\AqVg{s,-)\\Loodsln{e+ \\Vg{s,-)\\(^fids) 
q Jt Jt 



3 Proof of Theorem 1.1 



Since local existence of smooth solutions has been derived by N. Masmoudi [23] , here 
we only focus on the a priori estimate which is sufficient for proving Theorem ll.li As 
explained in [22] or [2^, to get the global existence we just need to control the L°°- 
norm of Vm, i.e., ||Vu||loo. Define the fiow associated with u by $(t,x), which means 
$ satisfies the ODEs, 

dt^{t,x) = u{t,^{t,x)), 

$(t = 0,a;) = X. 



(3.1) 



Step I Uniform estimates for tp and r with respect to t. 

Due to the special structure of the equation about ifjit^x^R), we have the follow- 
ing bounded estimates of if). For r > 1, multiplying the third equation of (11. ip by 



■000 


^o 


''1 


r 

B 



and integrating over B, then we get 



^ 



V^c 



dR + u-V / V^oo 
'b 



^ 



i^c 



dR 



4(r 



^c 



B 



Therefore, 



i'o 



i, 



-ipc 



{t,^{t,x),R)dR< / ^o 
Jb 



V'o 



-ipc 



Since the fiow is incompressible, then 

||V'IU?°,(£'-) < \\i'o\\L^{C-), 

To estimate r, we need a lemma. 



L-CLK/:--)) 



< 



'R 



i'oc 



{x,R)dR. 



MlLUcr)- 



dR. 



(3.2) 



Lemma 3.1. For any p such that pk > 1, it holds that 



1^1 



,bI-\Ri 

Proof. By the Holder's inequality 



-dR<C 



\i^\ 



p+i\ p+i 



B roo 



ijl 



£p+l. 



dR < C 

b1~\R\ - Jb{1- |i?|)i-^^/(^+i) ^Pj,(P^') 



<c 



Since pk > 1, the result follows. 



1 



B (1 - li^Di+i/p-'^ 






di? 



B i^oo 



D 



Noting that 



Tij{t, x)= {R^ ® VR^U)ip{t, X, R)dR, 
Jb 



one has 



|r(t,a;)| < C / |VijW| ■ |V(t,x, i?)|rfi? 
Jb 

^^^ I , ^^l.. -mt,x,R)\dR 
1 — \K\^ 



B 



f mt,x,R)\ 

-^Jb i-\R\ 

<c\m,x,-)\\cr, 



dR 



where in the last step we used Lemma [3TT1 
Hence, we have 

|kl|L°°(0,T;L2) < C||^o||l2(£'-), ||t||l°°(0,T;L°°) < C'UV^O ||Lg°(£'-)- 



(3.3) 



Step II A priori estimates for u. 

We need a useful lemma whose proof was established by Chemin and Masmoudi [7] 
(see also ilTI). 



Lemma 3.2 (Chemin-Masmoudi). Let v be a solution of the Navier-Stokes equations 
with initial data Vq G L^(M^), and an external force f G L^(0,T; C^^) fl L^(0,T; H~^): 



' dtv -Av + {v V)f + Vp = f, m M^ x (0, T) 
divf = 0, znM2x(0,T) 
v{t = 0, x) = Vq{x), in R^ 



(3.4) 



Then we have the following a priori estimates, 



and 



I'^IIli(0,T;C'1) 



+ (II^0||l2 + 



Il°°(0,T;L2) + 2||Vf ||^2(o^T;L2) < ||'yo||L2 + II /lli2(o,T;iJ-i)) ' 

< C Tsup ||A,t;o|U2(l - exp{-c22'^r}) 

+ sup / \\2-'^AJ{s)\\L^d^. 
q Jo J 



L'^{0,T;H 



-1 Vd 



L2(0,r;L2) 



Furthermore, if f & L^{0,T]C ^), then'ie > 0, there exists to(e) G (0,T) such that 



rllLi(*0,T';C'i) ^ <^- 



Particularly for our problem, since we have shown that r G L°°(0, T; L^)nL°°(0, T; L'^ 
applying Lemma 13. 2[ we know that 



u G L°°(0, T; L^) n ^^(0, T; ij^) n ^^(0, T; C^) 



and 



Ve > 0, 3 to(e) G (0, T), such that \\u\ 



Li(to,T;Cl) 



< e. 



(3.5) 



Step III Holder estimates for u 

For < t < T, choose some a satisfying < a < min {s — 2, 1}, define 



m{t,x) 



V^c 



B 



Agi;{t,x,R) 



-ipc 



dR= IIA, 



q'ywc 



{t,x) 



A{t) = sup ||m(s, ■)||(^i+c«, B{t) = sup ||r(s, ■)||(j,c« 

0<s<t 0<s<t 

D{t) = sup sup2"«||Arg(s,-)|U-. 
o<s<t <jez 

Here Ag is the frequency operator with respect to x. 

Before the detailed estimates, we will introduce an inequality for later use, which 
can be considered as an extension of Holder inequahty. 

Lemma 3.3. For any u E L^(R2) n C'1+"(R2), there holds that 

\\u (g) u\\^i+^ < C\\u\\l4 ■ ||m||^i+,, 
with some constant C independent of u. 
Proof. For any q E'L., using Bony's para-product decomposition [1], we have 

||Ag(M®M)|Uoo-2(5+°)'? <2 ^ ||A,(Viw®ApM)|Uoo.2(5+")'? 

|p-g|<5 
+ 5Z 5Z l|A5(Ap«®ArM)||Loo-2(5+")'? 
p>(j— 3 |p— r|<l 

= 2/i + h. 
Ii can be estimated as 

\p-q\<5 

<C J2 II Vi^lU^ ■ ||ApM||Loo ■ 2^^+"^" 

\p-q\<5 
|p-g|<5 

< C||m||/^4 ■ 1 1 M 1 1 (=,!+,,, 



here the first inequality is due to Lemma [2. II 



While I2 can be estimated as 

h = Y. Y. l|Ag(ApM®A,u)||Loo.2(^+")^ 

p>(j— 3 |p— r|<l 
p>ij— 3 |p— r|<l 

< C||m|||^i+„ ■ ||m||l4, 

here in the second inequahty we also used Lemma [2?n These above estimates complete 
the proof of Lemma 13.31 D 

First, applying Ag to the first equation of the FENE system, then we obtain 
dtAgU — vAAgU + VAgj9 = V ■ Ag(r — u^u), 

hence 

AgU = e^'^AgUo + I e^'^'-'^^ViAgV -{t-u® u))ds, 
Jo 

where P is the Helmoltz-Weyl projection operator. 

2^(i+°)||A,M|Uoo(t) 

< Ce-^"2''*||AgMo|U-2'?(^+") +C [ e-^'^2''(*-^)25(i+")||AgV ■ (r - u ® u)\\Lo^ds 

Jo 

< C||uo||ci+. +C f e-^22'(*-^) ■ 2-^(1+") (||V ■ A.rlUoo + ||V ■ Ag{u ® u)\\l^) ds 

Jo 

(3.6) 

Applying Lemma [2.11 we obtain 

'"*e-c2^'(*-)2'?(i+°)||V-A,r|Uoorfs 



'0 

< C re-^22'(*-^)2^(i+°)2'^||A,r|Uoocis 

< C I e-'=2''(*-^)22«-2"«||A„r||roods 



'0 

< CB{t). 



On the other hand, applying Lemma 12.11 again, we obtain 

Jo 
Jo 



< C I e-"22'(*-^)2i«||M®M|| .i^ds 





< C [ e 

Jo 



< C 



ic^+°' 



-C229(t-s)r,|g 



22'^||m||l4 • 1 1 M 1 1 (3,1+^ (is 

4 / rt 



I Il4 II ||4 7 

miLAMci+cds 



c22Ht-s)22q^, 



< c 



\u 



|2 IIV7 l|2 II ||4 r 



where we used Lemma 13.31 and interpolation inequality for the third and last inequality 
respectively. 

Taking the supreme of both sides of fl3.6p with respect to q, one gets 



By the Gronwall's inequality and the estimates in Step II, we have 

A{t) < C{\\uo\\ci+o. + B{t)) < C{1 + B{t)). 



(3.7) 



As a result of the relationship between r and ip and Lemma 13.11 

B{t) = sup sup2"''||Agr(s,•)||L- 
o<s<t q 



< sup sup 2""^ 

0<s« q 



\A,^\ ■ \VRU\dR 



B 



<C sup sup2-«||A,V^|U^^(^.)(s) 

0<s<t q 

= CD it) 
Combining the above inequality and (13. 7p . 

A{t)<C{l + D{t)). 



(31 



Step IV Holder Estimates of ip. 

The remaining part is to estimate the C'°-norm of ip. Take the operator Ag to the 
third equation, then 



dtAg^ + u ■ VAg^ = dwn{-W{u) ■ RA,^) + div^ (^00 V^ (^ 

+divK ([A„ -W{u)] ■ R^) + [m ■ V, A,]ij 



(3.9) 



10 



Multiply (El]) by r 



A,V 



r-2 



A,V 



^^ and integrate over B, 



dtK + u-VK 



4(r-l) 



B 



< f divR{[A,,-W{u)]-R^) 
Jb 









V^o 



r-2 



di? 






rdi? 



+ / |[wV,A,] 



A„V^ 



V'c 



r-1 



rdR 



= Ji + J2. 
By the Young's inequality and the Holder inequality, 



^0 



B 



R 



i'oc 



dR, 



Hence 



rr-2 I r^\\\„, V7 A l„/,ll ATr-1 



d,Nl + u ■ VNl < C\\ [A,, W{u)] ■ R^r^r ■ N;-' + C\\ [u ■ V, A^c^ ■ K 



g ' 



which implies that 

2''''^\\N,\\l^{t) < C2-^^' / \\[A„W{u)] ■ RtPWl^^cn ■ \\Ng\\l'-'i^)d^ 

Jo 



(3.10) 



+C2-''^' / ||b-V,AJ^|U: 



[.g°(£'-) ■ 



Ng\\l-J{s)dS. 



Note that 



\[A^,W{u)]-R^\ 



Then we get that 



h{y)[W{u){x) - W{u){x - 2~'^y)] ■ {Ri){x - 2"'y))dy 



<C2-^^W{u)\\a. / \h{y)\-\R^{x-2'^y)\dy, 



\[Ag,W{u)]-R^\\L^^cr)<C\\V 



u\ 



c° 



l-(z:-)-2""'<C||Vm||(.. -2- 



aq 



And by the Bony's para-product formula, 



|[u-V,A,]^| < J2 Yl l[ApM-V,A,]A,,^| 

P>9-3 |p-ij'|<l 

+ Yl \[Sq'-lU■V,Ag]Ag,^p\+ J2 |[A,'n-V,Ag]^g._i^| 

\q-q'\<5 \q-q'\<5 



11 



with each term having the following estimate: 
Since 

[ApM-V,A,]A,,V'(s,x,i?) 
= / h{y)[Apu{s,x)-Apu{s,x-2-''y)]-VAg,^{s,x-2-''y,R)dy, 

and 

VAg,ij{x) = [ 2'^^'h{2'i'y)Vi){x - y)dy = - f 2^'''Vh{2'''y)i;{x - y)dy 
then 

p>q—3 \p—g'\<l 

< C J2 Yl 2"^2-''||VApM|Uoo(s)-||VA,,^|Uc,(^.)(s) 

< CJ2 Yl 2-''2-''||VApM|Uoo(s)-2'''||^lUjo(^.)(s) 

p>q-3\p~q'\<l 

< C J2 2"1|VApiz|Uoo(s)-2"('?-P)||V^|Ugo(£.)(s) 

p>q-3 

< CA{s)-U\\L^^cn{s)<CA{s). 
Similarly, 

J2 2"%S,,_,u-V,Ag]Ag,ij\\L^^cr){s) 

\q-q'\<5 

< C Yl 2"«-'^||V5,,_in|Uoc(s)||VA,,^|Uoo(^.)(s) 

\g-q'\<5 

< C Y l|VM|Uo.(s)-2"'?'||A,,^|Ugo(£.)(s) 

\q-q'\<5 

< C\\Vu\\loo{s)-D{s). 

Y 2"''||[A,m-V,A,]5,,_i^|Uso(^.)(s) 
|g'-g|<5 

< C Y 2^'^-^VA,,u\\L^{s)-\\S,,.,Vij\\L^^cr){s) 

\q'-q\<5 

< CA{s)U\\L^^cn{s)<CA{s). 

Therefore, taking the supreme of fl3.10p with respect to q, by the Young's inequality 
and ([3JD, 

D{tY < f C[A{sY + D{sY]ds + C f \\Vu\\L^is)D{syds 
Jo^ Jo 

< [ C[l + D{sY] + C\\Vu\\L^D{sYds 
Jo 

12 



Then the Gronwall's inequahty imphes that 



Hence according to Lemma I2.2[ 

e + D{t) < C(D(0) + l)C,exp{C fli\\Vu\\Loo + l)ds} 

< CC,{D{0) + 1) exp{C(l + Jl \\u\\L2ds + eln(e + (t - QA{t)))} 

< CC4D{0) + 1) exp{C[l + eln(e + D{t))]} 
<CaiD{0) + l){e + D{t)f% 

where C* is some positive constant depending on the solution u on [0,t*]. Choosing 
e= J- 

D{T) < [Ca(D(0) + l)f. 
Then by fl3.8p . A{T) is bounded, which imphes that ||Vu||loo is bounded on [0,T] since 

II VuIIloo < C (11^11^2 + ||m||(^i+q) . 
4 Proof of Theorem 1.2 

The proof of Theorem 11.21 is similar to that of Theorem 11.11 we just give the sketch. 
As above, to get the global existence we only need to control ||Vm||loo, see [TO] . 

Step I Uniform estimates for ^ and r 

Integrate the third equation of (II. 3p on M, then 



dt I ijj{t,x,m)dm + u -W / ijj(t,x,7n)dm = 0. 

J M J M 

By the maximum principle for evolutionary equation, if) is always nonnegative. Hence 

\\^\LlnL^(L^M)){t) = ||^o|UinL-(Li(Af)) (4.1) 

lkl|L-'(0,T;L2) < C [\\lpo\\li(Li(^M)) + II^o|U2(l1(m)) j , (4.2) 

lkl|L°°(0,r;L-) < C{\\lpo\\l^(^L^M)) + 1) (4-3) 



Since s > | + 1 and M is a smooth compact manifold without boundary, 

Ml?^jh--(m)) < C||^o||l-(li(m))- (4.4) 

Step II A priori estimates for u 

13 



Applying Lemma \^72\ and the estimates (14.21) (14.31) . we get that 
and V e > 0, there exists to(e) ^ (0)^); such that ||M||i;i(-jQ t-c^) — ^■ 



Step III Holder estimates for u 

Denote 



J M 



ip(t,x,m)\'^d'm, 



A{t) = sup ||u(s, ■)||(3,i+^, B{t) = sup ||r(s, ■)||^„, 

0<s<t 0<s<t 

L>(t) = sup sup2"5||A/'g(s,-)IU-. 
o<s<t (jez 



As in section 3, we have the estimates 



A{t) < C(||mo||ci+. +5(t)) < C(l + 5(t)), 
and by g2D, Vg G Z, 

||A^g(t, OIU- < C'll^^li^, ■, ■)\\l^(lhm)) < <^||^o||l-(L1(M)), 



< 2 



-lz:/-l.,(2), 



/a/ /a/ ^mlHmlrij \{Hm,i^{s, X, mi) H m^^J {s , X, m2)) dmidm2 
+2'"i\^j^H-^-i'ff{m)H/\qip{s,x,m)dm 

<2"« y- I j H;^\H;^'j;^\mi,m2)Sp-iHi;{mi)\Hij{m2)dmidm2 
+2"' y, [ [ H;;,lH;;^l^if{mi,m2)ApH^{mi)Sp.iH^{m2)drmdm2 

+C2-5||Ar,(s,.)||^<3o 
< CD{s)\\H^\\l^^l2^m)){s) + CD{s) < CD{s) 

which implies that 

Bit) < CD{t). 

Therefore, 

A{t)<C{l + D{t)). 

Step IV Holder estimates for ip 



L°° 



L°° 



(4.5) 
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Take the operator H and A^ to the third equation, multiply by AgHip and integrate 
over M, then 

l-dt f \AgHiP\^dm + lu-V f \AgHiP\'^dm+ f \VgAgHiP\'^dm 

= / [u-V,Aq\Hilj-AqHipdm- / AqHdiVg {G{u,iIj)iIj) ■ AqHipdm 
Jm Jm 

= I [u-V, AglHip ■ AqH^dm - djUi [ HdiVg {ci^Agip) ■ HAgipdm 
Jm Jm 

+ f [djUi,Ag]H{c'iij) ■ VgAgHijdm+ I AgiVgUH^) ■ VgAgH^jdm 
Jm Jm 

+ [ Aq[HVgU,H-^]Hi)-VgAgH^dm 
Jm 

By the Young's inequality, we have 

< 2 / [u-V, Ag]Hip ■ AqHi/jdm 

Jo Jm 



ds 



+ 







djUi / HdiVg {d^Agilj)AgHipdm 



M 



ds 



+- / \\[djUi,Ag]H{d^ij)\\l^^^,^^^j^^ds 

^ U 

+ ? I \\Ag{VgUE^)\\ 







L^iL^M)) 



ds + - I \\Ag{\HygU,H-^]H^)\ 
4 Jo 



L-(L2(M)) 



/ ( Jl + J2 + ^3 + ^4 + ^5)C^S. 

Jo 



ds 



(4.6) 



As in section 3, applying fl4.5p . for every g G Z, 

22"''Ji(s) < CA{s)\\Hi,\\L^^^L2^^M)) ■ D{s) + C\\Vu\\loo ■ D{sf 
<C(D(s)2 + l)(||Vn|Uoc(s) + l). 

J2, J4 and J5 are estimated as in [9], 



< 



< 



M 



M 



HdiVg{d^Agip) ■ AgHipdm 
dwg {dlHAgi)) ■ AgH^dm 



M 



{[HdiYgd^, H-']H Ag^) ■ AgHijdm 



-(dwgC^)\AgHij\^dm 



M 



{[HdiVgcH, H-^]H Ag^) ■ AgHijdm 



< CWAgH^W 



L^(M)^ 
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which deduces that 

22^^2(5) < C22"5||Vm|Uoo(s) ■ N^{s) < C\\Vu\\loo{s) ■ D{sf. 

Using Bony's decomposition, 

2"^Ag{V gUHij)\\L^(^LHM)) 
<2-^ J2 ||ViV,W-A,i7^||^^(^,(,,)) 

|p-g|<5 

|p-g|<5 

p>q—3 |p— r|<l 
< C||VgW||L^(L2(M)) ■ SUp2°P||ApiJV^||L«=(L2(M)) 

p 

+ Csup2"P||ApVgW||L^(L2(jv/)) ■ ||i^V^||L-{L2{M)) 

p 
Combining the relationship between U and ip, 

2'^^.Ms) < C2^"1|i/^||iso(^2(M))(^) ■ N^is) < CD{sf. 
Similarly, 

p 
Since 

[d,u,, Ag]H{ci^ij) = [ h{y)[dMx) - dM^ - 2-~'iy)]H{c^ij){x - 2-''y, m)dy, 






'{L\M)) 



< CA{s)\ 
Therefore, taking the supreme of f l4.6p x 2"^°"^ with respect to t and q, 

D{tf <C f (||Vu||loo + 1)(1 + D{s)^)ds. 
Jo 

The remaining proof is just the same as that in section 3. We omit the details. 
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